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Abstract 

The effect of pressure and negative pressure on tlie modes of oscillation of a 
multi-electron bubble in liquid helium is calculated. Already at low pressures 
of the order of 10-100 mbar, these effects are found to significantly modify the 
frequencies of oscillation of the bubble. Stabilization of the bubble is shown to 
occur in the presence of a small negative pressure, which expands the bubble 
radius. Above a threshold negative pressure, the bubble is unstable. 
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Multielectron bubbles (MEBs) in liquid helium are fascinating entities expected to dis- 
play novel resonant behavior, the possibility of superconductivity, and sufficient electron 
density for Wigner crystallization and quantum melting. MEBs are bubbles inside the he- 
lium, containing only electrons that form a curved two-dimensional electron gas (2DEG) on 
the spherical surface of the bubble, the width of the thin spherical shell that conforms to 
the helium surface being of order 5-20 A compared to 70-80 A on a flat surface [§. 

New efforts to trap and localize MEBs for long periods of time @] have led to further con- 
sideration of the long term stability of bubbles. In this letter we discuss the use of pressure, 
and in particular negative pressure to stabilize MEBs against dynamic instabilities. 

The radius of the bubble depends on the enclosed charge. In a simplified model, valid 
for large bubbles, the radius is given by Rc = [e^ A^^/(167r ere)] where e is the electron 
charge, is the number of electrons in the bubble, a is the surface tension of helium and 
e is the dielectric constant of helium A single electron bubble has a radius of 17.2 
A while, for example, a bubble of 10^ electrons has a theoretical radius of 1 micron. 
There is some question about the static stability of an MEB: since the energy of the bubble, 
defined as the sum of the electrostatic and the surface tension energy, is proportional to 
N^'^ [g, clearly the energy of two bubbles with A^/2 electrons is lower than that of a single 
bubble with A^ electrons. Evidently fissioning is hindered by a formation barrier, since 
MEBs have been observed 0. Gravitational fields may fiatten very large bubbles and lead 
to instability [|l|. Dexter and Fowler showed that two-electron bubbles are unstable. 
Salomaa and Williams considered the dynamic stability against fissioning off of single 
electrons from large bubbles and found stability against this decay mode for bubbles with A^ 
greater than 15-20. It is straightforward to show that a positive pressure radially stabilizes 
a bubble, although angular modes can be unstable. With increasing negative pressure the 
bubble is first absolutely stabilized and then explodes, as we shall describe. Salomaa and 
Williams also considered the dynamic instability due to one of the surface oscillation modes, 
or ripplons, being soft (zero frequency) and found that this mode may be stabilized by 
the anharmonicity in the bubble's radial oscillation that results in a radius larger than Rq. 
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In this communication we study pressure related effects on the frequency of the modes of 
oscillation, on the equilibrium bubble radius, and on the stability of MEBs; we show the 
counterintuitive result that positive pressures can destabilize all higher angular modes, while 
negative pressures have a window of stabilization. We neglect gravity so that the MEBs are 
spherical. 

The frequencies of the modes of oscillation of a charged droplet were first calculated by 
Rayleigh P], and in the case of a charged bubble by Plesset and Prosperetti 0. We first 
set up a Lagrangian formalism to calculate the spherical ripplon modes. We then consider 
the effect of pressure on the static and dynamic properties of the bubble. The surface of 
the bubble is described by a function R{6, ip) that gives the distance of the surface from 
the geometrical center of the bubble, in the direction given by the two spherical angles 6, ip 
. This function can be written as R{9,(p) = _Rb + u{9 , ip) , where i?b is the angle-averaged 
radius of the bubble, and u{9, ip) describes the deformation of the surface from a sphere. 
This deformation can be expanded in a series of spherical harmonic deformations l^ml^, V^) 
with amplitude Q^^ 

oo e 

U{e, V) = Y,Y1 Q^rnYemie, ^) . (l) 

e=i m=-e 

In what follows, we shall assume that the amplitude of deformation is small, so that for all 
{£,m}, V£(£+l)Q,„//?b< 1. 

The kinetic energy T associated with the motion of the liquid helium surface can be 
derived from the velocity potential ip. The kinetic energy for incompressible flow of an 
inviscid fluid is given by 

T = ^jdipj d9,p{R^,e,ip) [n- V^(r,^,^)|,^J RlsmO, (2) 



with the density of helium p = 145 kg/m^, and n is the unit vector normal to the bubble 
surface. The coefficients of expansion in spherical harmonics of the velocity potential are 
then expressed as a function of the deformation amphtude Qim, leading to 
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The potential energy of the deformed bubble results from the surface tension, the pressure 
exerted on the bubble, and the electrostatic forces between the electrons of the spherical 
2DEG. The surface tension energy can be written as crS where o = 3.6x10^^ J/m^ is the 
surface tension of helium at zero pressure, and S is the area of the deformed surface : 

^ oo £ 

S = 4vri?^ + -^ ^(£2 + £ + 2)|g,„|^ (4) 

e=i m=-l 

Deforming the bubble will change its volume and perform work against the external 
pressure p of the helium liquid. The volume of the deformed bubble up to second order in 
the deformation amplitudes is given by §]: 

A OO I 

V=y/2^ + /?b$^ 5^ \Q£mt (5) 
1=1 m=-i 

The electrostatic energy of the electrons in the spherical 2DEG in the MEB can be derived 
by taking into consideration that the electrons are strongly confined in the direction per- 
pendicular to the helium surface (the binding energy is of the order of 10 K) and anchored 
to that surface, but free to move in the directions parallel to the helium surface. This leads 
to the following expression for the Coulomb part Uq of the potential energy M: 



JST2 2 JST2 2 ^ 



with e = 1.0572 the dielectric constant of helium. This expression is valid up to second order 
in the deformation amplitude Q^m and does not include exchange or correlation energies of 
the electron gas. 

Collecting the previous expressions (0,||,^,|D for the different energy contributions leads 
to the following expression for the Lagrangian iZbubbie = T — aS—pV—Uc '■ 



/:bubbie = - ^iraRl - —pRl - ^ 
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\Ql,n\^ ■ (7) 



The novel part in this Lagrangian as compared to previous treatments hes in the 

terms related to the pressure. We shall return to the harmonic solutions of this Lagrangian 
after discussing some static properties. 

The equilibrium radius of large bubbles, in the absence of deformations, Q^m = 0, can 
be found by minimizing the potential energy U = aS+pV+Uc as a function of the bubble 
radius. The potential energy is: 

C/W = -pii? + 4..flJ + ^j^^^ + ^ -0,3176-^ (8) 

For large bubbles, with ^ 10^, only the first three terms in the right-hand side of (H) 
play a role - these terms are also present in the Lagrangian (|^). The fourth term {m^ is the 
electron mass) is due to Shikin [|l| who makes the argument that there is a finite thickness 
(i <^ -R of the electron layer that should be taken into account, along with i?b, as a variational 
parameter to minimize the potential energy. The density functional calculations of Shung 



and Lin show that the exchange energy of the electron gas also plays a relevant role for 
small bubbles (N<1000), and that for practical purposes it is well approximated by adding 
an exchange term, the fifth term in the right-hand side of (|^). 

The equilibrium radius i?eq of the MEB is found as a function of the number of electrons 
and the exerted pressure, by minimizing the potential energy U{Rh) with respect to d and 
i?b- The value of -Rb which minimizes the potential energy ?7(i?b) is the equilibrium radius 
i?eq shown in Fig. 1 as a function of pressure for different numbers of electrons. The pressure 
reduces the equilibrium radius by a significant factor as compared to the zero pressure radius. 
Most of the change in i?eq occurs at low pressures. The behavior of Req as a function of 
pressure shown in Fig. 1 is similar for numbers of electrons up to 10* and larger. The graph 



extends to negative pressure |Tl[] increasing i?eq, since superfiuid liquid helium can sustain 
a substantial negative pressure. For any number of electrons in the bubble, there exists a 
critical negative pressure making the bubble unstable against runaway expansion. We found 
that the largest equilibrium radius that can be achieved before the critical underpressure is 
reached is approximately L5 times the equilibrium bubble radius at zero pressure, for any 



number of electrons from a few hundred to 10^. When the critical underpressure is reached, 
the only equilibrium radius is R^q ■—>■ oo. The critical underpressure appears to be inversely 
proportional to the number of electrons in the bubble, but this could not be shown explicitly. 

At fixed -Rb (the equilibrium radius), the part of the Lagrangian pertaining to the 
spherical ripplon modes represents a collection of harmonic oscillators in the coordinates 
Qim with ripplon frequencies 
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a{& + i + 2) + 2pR^-——^ 



(9) 



pR 

These frequencies are independent of the azimuthal index m. Fig. 2 shows the pressure 
dependence of the frequency of the spherical ripplon modes i = 1, 10. Note first that the 
£=1,2 modes have vanishing frequencies at pressures p ^ 0. This was discussed for zero 
pressure by Salomaa and Williams 0, who concluded that the £ = 2 mode can be dynami- 
cally stabilized if the effective radius of the bubble satisfies -Rb > -Rcq (the £ = 1 corresponds 
to uniform translation). However, the present treatment shows that this argument is no 
longer valid at increasing pressure. In fact, subsequent modes become unstable, as can be 
seen from Fig. 2 and expression (^. If the exchange and confinement energy terms in (H) 
are neglected (a reasonable assumption for bubbles with > 10^), the equilibrium radius 
satisfies 2pi?eq + 4(T = e'^N^ /ineRl , so that in this case 



P - 1 

i?b = i?eq ^UJe= \I—B^ W - 2) - 2pi?eq]. (lO) 



The pressure at which a mode £ > 2 becomes unstable is p = a{£ — 2)/ (2_Req). For = 10"^, 
a/Req is of the order of 100 mbar. Larger bubbles have even smaller critical pressures. The 
vanishing of the frequencies of these modes indicates that an instability occurs and that the 
quadratic approximation for the deformation is no longer valid. Note however, that a small 
negative pressure tends to stabilize the £ = 2 mode - this is an alternative to Salomaa and 
Williams' proposal for the stability of the £ = 2 mode, based on the assumption that the 
effective radius i?b > -Req- 

In this letter, we have shown that small negative pressures can stabilize a bubble against 
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dynamic instability, while positive pressures can drive all ripplon modes unstable. Since both 
negative and positive pressures are easily achievable experimentally, it will be interesting if 
bubbles can be created in a "stable" configuration that can be visually observed to study 
these predictions. In the above considerations we have considered the electrons to be a 
2DEG, ignoring Wigner crystallization. 
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FIGURE CAPTIONS 

FIG. 1. The equilibrium radius i? of a multi-electron bubble is shown as a function of 
pressure for bubbles for several values of A^. In the inset, the thickness d of the spherical 
electron shell in the bubble is shown as a function of pressure for the same numbers of 
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electrons. Both the equihbrium radius and the thickness were obtained by minimizing 
the potential energy (H). Negative pressure can be applied up to a critical underpressure, 
expanding the bubble. For = 3000 this critical underpressure is indicated with an arrow, 
as well as the maximum radius of a 3000 electron bubble. 

FIG. 2. The frequency of ripplon modes i = 1, .., 10 is shown as a function of external 
applied pressure (not vapor pressure) for = 1000. The leftmost points of the graphs start 
at the critical underpressure: for a pressure more negative than this value, the bubble is 
unstable against isotropic expansion. As the pressure is increased, more and more modes 
obtain a vanishing frequency. On the left the deformations (exagerated) are shown for a few 
of the modes. 



9 



This figure "MEBripfigl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0110569vl 



This figure "MEBripfig2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0110569vl 



